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Abstract 

Freely decaying two-dimensional Navier-Stokes turbulence is studied. The conservation of vor- 
ticity by advective nonlinearities renders a class of Casimirs that decays under viscous effects. A 
rigorous constraint on the palinstrophy production by nonlinear transfer is derived, and an up- 
per bound for the enstrophy dissipation is obtained. This bound depends only on the decaying 
Casimirs, thus allowing the enstrophy dissipation to be bounded from above in terms of initial data 
of the flows. An upper bound for the enstrophy dissipation wavenumber is derived and the new 
result is compared with the classical dissipation wavenumber. 

PACS numbers: 47.27.Gs, 47.27.Eq 
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In 1969, Batchelor 1 adapted Kolmogorov's equilibrium theory for three-dimensional (3D) 
turbulence to two-dimensional (2D) turbulence, on the basis of a phenomenologically anal- 
ogous property between the two systems. For a 3D fluid, decrease of the viscosity v is 
accompanied by increase of the mean-square vorticity, a consequence of the magnification of 
the vorticity by stretching of vortex lines, so that in the inviscid limit the energy dissipation 
is nonzero. For a 2D fluid, decrease of the viscosity enhances convective mixing, in which 
isovorticity lines get extended and brought closer to one another, 1 giving rise to increase of 
the mean-square vorticity gradients (twice the palinstrophy) , so that in the inviscid limit the 
rate of enstrophy (half the mean-square vorticity) dissipation can approach a finite value \. 
On the basis of this analogy, Batchelor 1 applies the familiar arguments of the Kolmogorov 
equilibrium theory for the small-scale components of 3D turbulence to the 2D case, where 
the roles of the energy and energy dissipation in the original theory are played by the enstro- 
phy and enstrophy dissipation. This means that the statistical properties of the small-scale 
components of the turbulence depend only on the two dimensional parameters x an d v. The 
enstrophy dissipation \ is thus an important dynamical quantity in Batchelor's theory. One 
of its prominent role is in the expression of the enstrophy spectrum Z(k) of the so-called 
enstrophy inertial range, which is presumably formed when an initial enstrophy reservoir 
spreads out in a virtually inviscid region of wavenumber space: 

z(k) = c x 2/3 k-\ (1) 

where C is a universal constant and k is the wavenumber. Another important role of \ is 
in the determination of the dissipation wavenumber k v : 

X 1/6 

which presumably marks the end of the enstrophy inertial range, around which the enstrophy 
is most strongly dissipated. In both (JTJ) and (J2J), X is a finite (but otherwise undetermined) 
parameter. 

It is desirable to have a quantitative knowledge of x (and hence of k u ), not only for 
its role in the Batchelor theory but also for further analyses of the turbulence, beyond 
the usual dimensional arguments. Even for finite Reynolds numbers, the determination 
of x is highly non-trivial. In the limit u — > 0, this problem does not seem to become 
more tractable. Ideally, x can be determined if the extent to which the production of 
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palinstrophy by convective mixing is fully understood. This letter takes a direct approach 
to this problem by deriving a rigorous upper bound for the nonlinear term representing the 
palinstrophy production rate. Equating this bound to the viscous dissipation term yields a 
constraint, from which upper bounds for \ and for the enstrophy dissipation wavenumber k& 
(to be defined later in this letter) can be derived. These bounds are found to be completely 
described in terms of initial data of the flows. The derived enstrophy dissipation wavenumber 
is consistent with the classical prediction of k v given by (J2J), in the sense that they both have 
the same functional dependence on v. A novel result of this study is that x an d k^ can be 
estimated in terms of initial data of the turbulence, while they (more precisely x an d k v ) 
are essentially undetermined in the Batchelor theory. 

In the vorticity formulation, the freely evolving 2D Navier-Stokes equations governing 
the motion of an incompressible fluid confined to a doubly periodic domain are 



where £(x, t) is the vorticity, J(6, $) = B x 'dy — 9y'& x , v the kinematic viscosity, and ip(x, t) the 
stream function. The vorticity is defined in terms of the stream function and of the velocity 
v by £ = A^ = n ■ V x v, where n is the normal vector to the fluid domain. Equivalently, 
v can be recovered from ip and £ by v = (— ip y , ip x ) = (— A -1 ^, A -1 ^). 

An importance property of the advective nonlinear term in (JHJ) is that it conserves the 
kinetic energy and an infinite class of integrated quantities, known as Casimirs, including 
the enstrophy. The latter conservation law is attributed to the fact that vorticity is conser- 
vatively redistributed in physical space by the advective transfer. While the conservation of 
energy and enstrophy imposes strict constraints on turbulent flows and has been explored 
in the literature to a great extent, 2-19 the conservation of Casimirs, other than the enstro- 
phy, seems to render little additional knowledge of the flows and has received much less 
attention. 20-22 In the present case of unforced dynamics, it is well known that a wide class 
of these Casimirs decays under the action of viscosity. Here our main interest is in the 
following Casimirs: (|£| p ), where (•) denotes a spatial average. By taking the time derivative 
of (|£| p ) and using (J3J) one obtains 



d t £ + J(^,£) = ^A£, 



(3) 



d 
dt 



v(\tr 2 £J(M) + 'v(\t\ p - 2 tA£) 

-^(p-l)(|£r 2 |V£| 2 ), 



(4) 
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where the second equation is obtained by integration by parts, upon which the nonlinear term 
identically vanishes. For p — 2, Eq. (@J) governs the decay of (|£| 2 ) (twice the enstrophy), for 
which the dissipation term becomes 2%, which is the subject of this study The right-hand 
side of (j3J) is negative for p > 1. Hence decays in time, for p > 1. It follows that 

(\m\ p ) < (ie(o)H, (5) 

for p > 1 and t > 0. In particular, in the limiting case p — > oo, one has 

iie(*)L<ne(o)L, (6) 

for t > 0, where ||£|loo denotes the L°° norm of £. 

Now the main result of this letter can be readily derived. By multiplying Q by A£ and 
taking the spatial average of the resulting equation, one obtains the equation governing the 
evolution of the palinstrophy (|V£| 2 )/2: 

lj t m\ 2 } = (^j^,o)-vm\ 2 } 

< (leiawgiv&i + ivv v iive„i)) - Hm 2 ) 

< iieiL (iv^r + iv^i 2 ) i/2 (ive,i 2 + ive/> i/2 - urn?) 
<neiL(ier) i/2 (iAei 2 ) i/2 -^(iAer). (?) 

In (J7J), the second equation is obtained via the two elementary identities 

A Jty, = AO + 2 J(^, + 2 J(^, £,) (8) 

and 

(A^,0) = -(^,AO). (9) 

The Holder inequality is used in the fourth step, and the last step can be seen by expressing ip 
(and £) in terms of Fourier series. The triple-product term (|0 2 ) 1 ^ 2 (| A^ 2 ) 1 ^ 2 represents 
an upper bound for the palinstrophy production rate. It can be seen that the palinstrophy 
necessarily ceases to grow as its dissipation u(\ A£| 2 ) reaches this bound. It follows that as the 
palinstrophy grows to and reaches a maximum (d(|V£| 2 ) /dt > 0), the following inequality 
necessarily holds 

(|A£| 2 ) 1/2 < ^^ V2 , (10) 
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Since (|V£| 2 ) can be bounded from above in terms of (|A£| 2 ), ineq. (fTUjl can be used to 
derive an explicit upper bound for the palinstrophy. By Holder inequality one has 18 

(m 2 )>^§f, (ii) 

where the inequality sign ">" can become "3>" (see below). Substituting (fTTJ) into (fTUj) 
yields 

(iv£i 2 > < SoM!l. (12 ) 

It follows that the enstrophy dissipation x is bounded from above by 

x < iieiioo <i^r>- as) 

It is notable that the upper bound for \ in (|T3|) is expressible in terms of two decaying 
Casimirs, namely and (|£| 2 ), so that it can be bounded from above in terms of initial 
data of the flows. More accurately, x can t> e bounded from above in terms of the initial 
vorticity field only. In passing, it is worth mentioning that since H^H^ and (|£| 2 ) are intensive 
quantities, i.e. independent of the domain size, the constraint (fT3|) is size-independent. 

Let us denote by k d and k D the wavenumbers defined by (| V^| 2 ) 1 / 2 /(|^| 2 ) 1 / 2 and 
(| A^| 2 ) 1 / 2 /(| V£| 2 ) 1//2 , respectively. For "regular" spectra, k d [k D ) specifies where, in 
wavenumber space, (|V£| 2 ) ((|A£| 2 )) is mainly distributed. In other words, kd (kn) spec- 
ifies where, in wavenumber space, the enstrophy (palinstrophy) dissipation mainly occurs. 
By "regular" it is meant that the enstrophy is not highly concentrated in any particular 
regions of wavenumber space that would result in severe steps in the enstrophy spectrum. 
By these dissipation wavenumbers satisfy k d < ko, and the inequality sign "<" can 
become "<C". This can be realized if the palinstrophy spectrum around kd does not fall 
off so steeply. More quantitatively, by the definition of kd, one can expect the palinstrophy 
spectrum around k d to be shallower than fc -1 . (Because, otherwise, most of the contribution 
to (|V£| 2 ) would come from k < k d , making the ratio (| V^| 2 ) 1 / 2 /(|^| 2 ) 1 ^ 2 significantly lower 
than kd, a contradiction to the very definition of kd-) This means that the spectrum of 
(| A£| 2 ) around kd is shallower than k 1 . Hence, most of the contribution to (|A£| 2 ) can come 
from k ^> kd if the palinstrophy spectrum beyond kd becomes steeper than k~ x and falls off to 
g rac [ ua iiy. This allows for the possibility ko ^> k d to be realized. In any case, k d should 
be well beyond the end of the enstrophy range, i.e. J k>kd Z{k) dk/ J k<kd Z{k) dk w 0, and kr> 



should be well beyond the palinstrophy range, i.e. f k>kD P(k) dk/ f k<kD P(k) dk ps 0, where 
Z(k) and P(k) are the enstrophy and palinstrophy spectra, respectively. In other words, the 
enstrophy spectrum around kd should be steeper than fc -1 , and the palinstrophy (enstrophy) 
spectrum around kr> should be steeper than k~ x (k~ 3 ). 

Our primary concern is an estimate of kd when (|V£| 2 ) achieves a maximum. It is likely 
that kd achieves a global maximum then. Equation (JJJ) and the subsequent equations (fTUj) 
and (JT2jl imply 

k d k D < Si (14) 
v 



fcd < ( imioo i < / Kik ) (15) 



and 

, 1/a < pie(Q) 11 Nl/2 

It is interesting to compare the present result with the classical dissipation wavenumber 
k u . To this end, let us express \ — u (\^i\ 2 ) m the form \ — u ^d{\i\ 2 )i so that (J2J) can be 



rewritten as 



x i/e fk d (\e) i/2 ^ 1/3 

z/ 1 



^ = ^=(^L^) . (16) 



It follows that 

ki mY 12 



(17) 



fc d v 

It may be assumed that H^H^ and (|£| 2 ) 1//2 are of the same order of magnitude. In fact, one 
can even have the exact equality HClloo = (|£| 2 )^ 2 f° r some simple cases. For example, for a 
vorticity field of a single Fourier mode, this equality trivially holds. In any case, both 
and (|£| 2 ) 1//2 decay under the action of viscosity, so that if they are comparable initially, then 
one can expect them to be comparable subsequently. Hence one can deduce from ((Hj) and 
(PZD that 

k 2 d k D « kl (18) 

It follows that kd < k u , and kd can be significantly lower than k v if kd^iku. 

As pointed out by a referee, the result (JT5|l can be derived from more physical consider- 
ations, based on the relation 

d t {i{l)i{2)) + 9 r (<K(l)£(2)> « vd 2 Ml)t{2)), (19) 



where 1 and 2 stand for x\ and x 2 , respectively, r = \x 2 — X\\, and 8v denotes the fluid speed 
increment across r. In the limit r — > , if one assumes £ ~ 5v/r then one can recover (JTSJl 
by balancing the nonlinear and viscous terms in (|19|). 

In conclusion, this letter has derived a rigorous constraint on the palinstrophy production 
rate by nonlinear transfer, from which upper bounds for the enstrophy dissipation x and for 
the enstrophy dissipation wavenumber kd have been deduced. These bounds are expressible 
in terms of the vorticity supremum the mean-square vorticity density (|£| 2 ), and the 

viscosity v. These quantities are "intensive" quantities, i.e. independent of the domain size, 
making the derived upper bounds size-independent. Moreover, these bounds are completely 
determined by the initial vorticity field and v since both ||£|| and (|£| 2 ) decay under the 
action of viscosity. The upper bound k d < ||£(t)||^ 2 /V 1 / 2 < ||£(0)||^ 2 jv x l 2 is consistent 
with the classical dissipation wavenumber k v = x 1 ^ 6 /^ 1 ^ i n the sense that they both have 
the same functional dependence on v. A novel result of this study is that both x an d kd 
are bounded in terms of the initial vorticity field and the viscosity whereas x an d k v are 
essentially undetermined by the classical theory. 
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